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^\ • Abstract 

| Coalescing compact binaries have been pointed out as the most promising source of grav- 

j^Qi itational waves for kilometer-size interferometers such as LIGO. Gravitational wave signals 

r-{ ' are extracted from the noise in the detectors by matched filtering. This technique performs 

, really well if an a priori theoretical knowledge of the signal is available. The information 

known about the possible sources is used to construct a model of the expected waveforms 
(templates). A common assumption made when constructing templates for coalescing com- 
pact binaries is that the companions move in a quasi-circular orbit. Some scenarios, however, 
predict the existence of eccentric binaries. We investigate the loss in signal-to-noise ratio 
induced by non-optimal filtering of eccentric signals. 

co ; 

q ; 1 Introduction 

°o ' 

Coalescing compact binaries have been pointed out as the most promising source of gravitational 
0^ , waves for the LIGO /VIRGO /TAMA/GEO interferometers [| §. These binaries typically have 

formed a long time ago, giving them time to radiate most of their eccentricity away. The 
templates (model of the radiation) needed for matched filtering are thus constructed according 
to this assumption. Gravitational waves emitted by a circular binary will be explicitly searched 
i—t , for in the output of the detectors, but not gravitational waves emitted by eccentric binaries. 

The scenario we have in mind allows for the formation of young eccentric binary systems, young 
, enough that they did not have had time to be fully circularized by the radiation reaction. For 

example, the collapse of a dense Newtonian globular cluster can lead to the formation of a copious 
number of eccentric binaries via two- and three- body encounters || ||. 

These eccentric binaries will emit strongly in the frequency band of the LIGO interferometers. 
It may seem that these eccentric binaries can be dealt with by incorporating adequate templates 
in the bank of templates already available, but this may prove to be inefficient. The addition of 
new templates has two undesirable effects: It adds to the already heavy computational burden 
associated with data processing and it increases the probability of false detection (mistaking 
the noise in the detector for a signal). A better solution might be to search for these eccentric 
signals with the circular templates, and once a signal is concluded to be present, to extract the 
information using eccentric templates. 

For this to be possible, the circular templates have to follow the phase of the eccentric signals 
very well. To assess the quality of the circular templates at modeling eccentric signals, special 
detection tools are needed. 
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2 Matched filtering as a detection method 

Gravitational wave signals are very weak and at best they will be of the same order of magnitude 
as the noise in the detectors. This motivates the general belief that matched filtering will be 
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needed to extract the signals from the noisy output of the detectors [[0. When the signals are of 
known shape, this technique produces the highest signal-to-noise ratio ]!|. Suppose a gravitational 
wave h(t) reaches the detector. The output of the detector o(t) is then a superposition of the 
useful signal h(t) and the noise n(t). In matched filtering, the signal is extracted by using a 
theoretical template (theoretical model) that mimics the signal as well as possible; we call this 
template m(t,Sl). The vector ft denotes the parameters that characterize the template. If the 
template were a perfect copy of the signal, the parameters would represent the real parameters 
of the source, such as its mass and distance from earth. If the templates are not a perfect 
approximation to the real signal, the parameters H represent phenomenological parameters. 

If instead of working in the time domain we work in frequency space, we can introduce 
the natural inner product of matched filtering. For two functions a(t) and b(t) with Fourier 
transforms 5(/) and b(f), the inner product is defined as|| 

m = 2 r d/ a -ww+yw> , (1) 

JO ^n(I) 

where a " * " denotes complex conjugation and S n (f) is the one-sided spectral density of the 
detector's noise. In terms of this inner product, the average signal to noise ratio is@] 

(P) = / {mim) ■ (2) 

In practice, the set of parameters ft is varied until a maximum of the signal-to-noise ratio 
is found. This maximum is the signal-to-noise ratio achievable by using rh(f, fi) as a template. 
The signal-to-noise ratio of equation (||) does not give any information about the quality of the 
templates or, equivalently, how well the template models the signal. The Schwartz inequality 
provides an answer to this questionjQ. The absolute maximum the signal-to noise ratio can take 
is achieved when the template is a perfect match of the signal, and the parameters f2 correspond 
to the parameters of the source (rh(t, ft) = h(f)). The optimal SNR is (7) 



(p) max = VW0- (3) 



By dividing the signal-to-noise ratio (equation (|2|)) with the value achieved by optimal filtering 
(equation (^)), we construct the ambiguity function A(H); 

Am = . ( 4) 

y/(m(n)\m(n))(h\h) 

This function takes values between and 1. It is equal to 1 when the optimal template is used. 

The value of the parameters $7 can be varied until A(St) is maximized. The maximum value 
of the ambiguity function is the fitting factor: 

FF = max A(fl) . (5) 

The fitting factor is a direct measure of the template's quality since it can be related to the loss 
of event rate, i.e. the number of events missed by using an inappropriate set of templates. This 
loss is calculated according to 1 — i^F 3 [||. For example, if the fitting factor is 0.8, then 48.8% 
of the events would be mistaken for noise. We adopt a threshold of FF = 0.9 for the present 
work. This corresponds to a loss in event rate of 27%. 
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3 The gravitational waveforms 



We calculate the waveforms for both circular and eccentric binaries in the quadrupole approxi- 
mation. In this approximation the waveforms are given by 13] 



2 d 2 



hf/ = —— I I, 
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(6) 



where R is the distance between the source and the observer, Itj is the source's quadrupole 
moment and the superscript TT reminds us that gravitational waves are traceless and live in the 
plane transverse to the direction of propagation. 

For eccentric binary systems, the waveforms areS 



x = - h yy = 2 (! + cos2 #o) cos2(</? - ip ) 

Hp 

(1 + cos 2 o ) cos(ip — 2ip ) + i cos(3(^ — 2<p ) ) + sin 2 9„ cos y - 



+ e 2 (l + cos 2 O ) cos2^ -I- sin 2 



(7) 
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cos 9 < 4 sin 2 (</? — 2(/? ) 



+ e [5sin(y; — 2ip a ) + sin(3^ — 2<p )} — 2e 2 sin2(^ > , 



(8) 



where £? and <^ are the two angles defining the location of the observer with respect to the 
orbital plane, fi is the reduced mass, and p and e are defined in terms of the turning points (r±) 
of the Newtonian orbit as 
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The eccentric waveforms oscillate at once, twice and thrice the orbital frequency, whereas the 
circular waveforms oscillate only at twice the orbital frequency. We parameterize our binaries 
by specifying the two masses, and the eccentricity e G they have when they first enter the LIGO 
frequency band (at 40 Hz). 

The Fourier transform of the waveforms is calculated numerically for the eccentric signal, 
and obtained through the stationary phase approximation for the circular templates ]ft0||. Once 
these Fourier transforms are known, it is straightforward to calculate the optimal SNR (equation 
), build the ambiguity function (equation (^)), and maximize it over the different parameters 
of the templates to get the fitting factor (equation (||)). 

The results for the signal-to noise ratio and the fitting factor are displayed in figure (jl]) . The 
signal-to-noise ratio for an eccentric signal is higher than the ratio obtained for an equivalent 
circular binary. This means that if both binaries are located at the same distance i?, the eccentric 
binary will emit stronger radiation and will be easier to detect if optimal filters are used. On the 
other hand, the fitting factor decreases as the initial eccentricity is increased. The circular tem- 
plates fail to model the eccentric signal properly. The good news is that circular templates are 
still accurate enough to detect some eccentric signals. For example, a neutron star binary system 
will be detected as long as its initial eccentricity does not exceed 0.13. If the total mass of the 
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Figure 1: Figure A: The ratio of the optimal signal-to-noise ratio (p(e D )) to the signal-to-noise 
ratio ((p(Q))). The figure shows that eccentric binary systems will be easier to detect if they are 
explicitly searched for at the output of the interferometers. 

Figure B: The fitting factor as a function of e Q for various binary systems. Two trends are 
apparent. The first one is the net decrease in the fitting factor as e G increases, while the total 
mass of the binary is held fixed. The second one is the increase in the detection probability 
when the total mass of the binary increases. The various binaries studied are labeled by the two 
masses of the companions; they are given in units of the solar mass. 



system is increased, the detection probability increases as well. For example, for a system of two 
8.OM0 black holes, the initial eccentricity can be as high as 0.33. This trend is explained in the 
following way As the total mass of the system increases, the radiation it emits is stronger and 
the system coalesces in a shorter time. The shorter the signal, the less opportunity the circular 
templates have to go out of phase with it. This is a good news, because the higher the mass is, 
the easier it is to detect the signal because it is stronger. Thus, the LIGO interferometers should 
be able to detect radiation from some eccentric binaries, those with large masses and relatively 
low eccentricities p"H . 
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